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Abstract 


The objective of this thesis is to study the existence and uniqueness of the solution of a nonlinear 
integro-differential equation of Volterra of second kind when the derivative of the solution 
appears under the sign of integration in a non-linear way, which has been defined without 
a singular kernel. Then, using Matlab and C ++ we get the approximate solution and error 
estimates. During this study, we deemed it necessary to study a problem of Fractional Boundary 
value in the presence of this derivative, it leads to the study of a linear integral Volterra- 
Fredholm equation of the second type and this by defining a new fractional integral as inverse 


of the conformable fractional derivative of Caputo. 
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Résumé 


L’objectif de cette these est d’étudier l’existence et l’unicité de la solution d’une équation 
intégro-différentielle non linéaire de Volterra de second type lorsque la dérivée de la solution 
apparait sous le signe de l’intégration de maniére non linéaire, qui a été définie sans noyau 
singulier. Ensuite, en utilisant Matlab et C ++, nous obtenons la solution approximative et les 
estimations d’erreur. Au cours de cette étude, nous avons jugé nécessaire d’étudier un probleme 
fractionnaire aux limites en présence de cette dérivée, ce dernier conduit a l’étude d’une équation 
Volterra-Fredholm intégrale linéaire du second type et ce en définissant une nouvelle intégrale 


fractionnaire comme inverse de la dérivée fractionnaire conformable de Caputo. 


Mots clés: 
Equation intégro-différentielle non linéaire de Volterra, Dérivée fractionnaire sans noyau sin- 
gulier, Convergence, Solution analytique, Existence et unicité de la solution, Problemes aux 


limites fractionnaire, fonction de Green, intégrale fractionnaire. 
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Introduction 


This thesis is the outcome of my research during my Ph.D. study at Guelma University. The 


principal materials in the thesis are based on the following articles from this period : 


Paper I: Moumen Bekkouche M., Guebbai H. & Kurulay M. Analytical and numerical 
study of a nonlinear Volterra integro-differential equations with conformable fractional 
derivation of Caputo. Annals of the University of Craiova - Mathematics and Computer 


Science Series(2020). https://doi.org/. 


Paper II: Moumen Bekkouche M., Guebbai H. & Kurulay M.On the solvability fractional 
of a boundary value problem with new fractional integral. J. Appl. Math. Comput. (2020). 
https://doi.org/10.1007/s12190-020-01368-x. 


Paper III: Moumen Bekkouche M., Guebbai H. Analytical and Numerical Study for 
an Fractional Boundary Value Problem with conformable fractional derivative of Caputo 
and its Fractional Integral. J. Appl. Math. Comput. Mech., JAMCM(2020). https: 
//doi.org/ 


Paper IV: Moumen Bekkouche M., Guebbai H. & Kurulay M., S. BenmahmoudA 
new fractional integral associated with the Caputo—Fabrizio fractional derivative. Ren- 
diconti del Circolo Matematico di Palermo Series 2 (2020). https://doi.org/10.1007/ 
$12215-020-00557-8. 


The organization of this thesis is as follows : The thesis begins in Chapter 1 that contains 
a brief history about the appearance of the concept of fractional calculus and a presentaion 
of some fractional derivative. In the sequel, we describe some basic concepts and give some 
information about direct and inverse problems. Afterward, we collect the results obtained from 


three articles above and present them respectively in Chapters 2, 3 and 4. 


The next three chapters are structured on fractional calculus as follows: We start in Chapter 
2, by studying an integro-differential equation of Volterra when the derivative of the solution 
appears under the sign of integration in a non-linear way. The derivation is conformable frac- 
tional of the Caputo type, we prove two theorems of existence and uniqueness for the solutions 
of this problem in a Banach space specific to the derivative of Caputo. The analytical study is 
followed by a complete numerical study. 

In Chapter 3, we introduce a new definition of fractional integral as an inverse of the 
conformable fractional derivative of Caputo. By using these definitions, we obtain the basic 
properties of those fractional derivative and its fractional integral. Finally, we solve a class of 
fractional boundary value problems as a theoretical application, and we use Matlab to solve this 
class of fractional boundary value problems. In Chapter 4, we consider a class of a fractional 
boundary value problem with conformable fractional derivation of the Caputo type. We obtain 
existence and uniqueness results for this problems based on the new definition of fractional 
integral. Therefore, the proofs are based upon the reduction of the problem to a equivalent 
linear Volterra-Fredholm integral equations of the second kind and this analytical study is 


followed also by a complete numerical study. 
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Introduction and Preliminaries of 


Fractional Calculus 


Summary 
1.1 History, Definitions and Applications .................. 11 
1.1.1 Griinwald—Letnikov derivative .................. 15 
1.1.2 Riemann-Liouville fractional derivative ............. 15 
(Re uk, a ee ee 16 
1.1.4 A new fractional derivative of Caputo. ............. 16 
1.2 Differentiation under the Integral Sign ................. 21 


For the convenience of the reader, in this chapter we present here the necessary definitions and 
fundamental facts of the fractional calculus theory in this chapter. These definitions can and 


fundamental facts be found in the recent literature. 


1.1 Fractional Calculus: History, Definitions and Appli- 
cations 


This part is aimed at the reader who wishes to learn about Fractional Calculus and its possible 
applications in his/her field(s) of study. The intent is to first expose the reader to the concepts, 


applicable definitions, and execution of fractional calculus (including a discussion of notation, 


iBi 


CHAPTER 1. FRACTIONAL CALCULUS 


operators, and fractional order differential equations), and second to show how these may be 
used to solve several modern problems. 

The traditional integral and derivative are, to say the least, a staple for the technology 
professional, essential as a means of understanding and working with natural and artificial 
systems. Fractional Calculus is a field of mathematics study that grows out of the traditional 
definitions of the calculus integral and derivative operators in much the same way fractional 
exponents is an outgrowth of exponents with integer value. Consider the physical meaning of 
the exponent. According to our primary school teachers exponents provide a short notation 
for what is essentially a repeated multiplication of a numerical value. This concept in itself 
is easy to grasp and straight forward. However, this physical definition can clearly become 
confused when considering exponents of non integer value. While almost anyone can verify 
that «2 = x-x-2x, how might one describe the physical meaning of x*4, or moreover the 
transcendental exponent x!/4. One cannot conceive what it might be like to multiply a number 
or quantity by itself 3.4 times, or : times, and yet these expressions have a definite value for 
any value «x, verifiable by infinite series expansion, or more practically, by calculator. Now, 
in the same way consider the integral and derivative. Although they are indeed concepts of a 
higher complexity by nature, it is still fairly easy to physically represent their meaning. Once 
mastered, the idea of completing numerous of these operations, integrations or differentiations 
follows naturally. Given the satisfaction of a very few restrictions (e.g. function continuity), 
completing n integrations can become as methodical as multiplication. But the curious mind 
can not be restrained from asking the question what if n were not restricted to an integer 
value? Again, at first glance, the physical meaning can become convoluted (pun intended), but 
as this report will show, fractional calculus flows quite naturally from our traditional definitions. 
And just as fractional exponents such as the square root may find their way into innumerable 
equations and applications, it will become apparent that integrations of order 1/2 and beyond 
can find practical use in many modern problems. 

The concept of fractional calculus is popularly believed to have stemmed from a question 
raised in the year 1695 by Marquis de L’Hopital (1661-1704) to Gottfried Wilhelm Leibniz 
(1646-1716), which sought the meaning of Leibniz’s (currently popular) notation ou for the 
derivative of order n € No := {0,1,2,3,---} when n = 1/2 (What if n = 1/2? ). In his 


” 


reply, dated 30 September 1695, Leibniz wrote to L’Hopital as follows: ”... This is an apparent 


paradox from which, one day, useful consequences will be drawn. ...” 
Subsequent mention of fractional derivatives was made, in some context or the other, by 


(for example) Euler in 1730, Lagrange in 1772, Laplace in 1812, Lacroix in 1819, Fourier in 
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1822, Liouville in 1832, Riemann in 1847, Greer in 1859, Holmgren in 1865, Griinwald in 1867, 
Letnikov in 1868, Sonin in 1869, Laurent in 1884, Nekrassov in 1888, Krug in 1890, and Weyl 
in 1917. In fact, in his 700-page textbook, entitled ”Traite du Calcul Differentiel et du Calcul 
Integral” (Second edition; Courcier, Paris, 1819), S. F. Lacroix devoted two pages (pp. 409-410) 


to fractional calculus, showing eventually that 


Leibniz’s response, based on studies over the intervening 300 years, has proven at least half 
right. It is clear that within the 20th century especially numerous applications and physical 
manifestations of fractional calculus have been found. However, these applications and the 
mathematical background surrounding fractional calculus are far from paradoxical. While the 
physical meaning is difficult (arguably impossible) to grasp, the definitions themselves are no 
more rigorous than those of their integer order counterparts. Understanding of definitions and 
use of fractional calculus will be made more clear by quickly discussing some necessary but 
relatively simple mathematical definitions that will arise in the study of these concepts. These 
are The Gamma Function, The Beta Function, The Laplace Transform, and the Mittag-Leffler 


Function. 


Definition 1.1 (The Gamma Function). the gamma function (represented by I. , the capital 
letter gamma from the Greek alphabet) is one commonly used extension of the factorial function 
to complex numbers. The gamma function is defined for all complex numbers except the non- 


positive integers. For any positive integer n, 


I(n) =(n—-1)!. 


Derived by Daniel Bernoulli, for complex numbers with a positive real part the gamma function 


is defined via a convergent improper integral: 


Ft) = ate "dz, forallzéR 
0 
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| 
¢ 


Figure 1.1: The gamma function along part of the real axis 


The gamma function then is defined as the analytic continuation of this integral function 
to a meromorphic function that is holomorphic in the whole complex plane except the non- 
positive integers, where the function has simple poles. The gamma function has no zeroes, so 


the reciprocal gamma function 1/T is an entire function. 


The first mention of a derivative of arbitrary order appears in a text. S. F. Lacroix [1819, 
pp. 409 - 410] devoted less than two pages of his 700-page text to this topic. He developed a 
mere mathematical exercise generalizing from a case of integer order. Strating with y = 2”, m 


a positive integer, Lacroix easily develops the n—th derivative 


Using Legendre’s symbol for the generalized factorial (the gamma function), he gets 


d”y T(m +1) 
dx” T(m—n-+1) 


m—-n 


He then gives the example for y = x and n = 5, and obtains 


1 
2? 


d'/*y = off 
dx? /q 
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It is interesting to note that the result obtained by Lacroix, in the manner typical of the 
classical formalists of this period, is the same as that yielded by the present-day Riemann- 
Liouville definition of a fractional derivative. Lacroix’s method offered no clue as to a possible 
application for a derivative of arbitrary order. Now we have several definitions of the fractional 


derivative, we present: 


1.1.1 Grtnwald—Letnikov derivative 


Griinwald—Letnikov derivative is a basic extension of the natural derivative to fractional one, 
It was introduced by Anton Karl Griinwald in 1867, and then by Aleksey Vasilievich Letnikov 


in 1868. Hence, it is written as 


Se Qa 
De f(t) = tim 1" || stm) 
m=0 m 


1.1.2 Riemann-Liouville fractional derivative 


Riemann-Liouville fractional derivative acquiring by Riemann in 1847 is defined as follows. 
gee d\" ( 5-(n-a) 
per) =(2) (ie) £0) 


a ical (5) [ “ae (n = [a] +1,t >a) 


where a@ > 0 this operator is an extension of Cauchy’s integral from the natural number to real 


one. Based on the fractionalization algorithm. In addition, according to the above relation, if 


0<a<1 then the Riemann-Liouville operator reduced to 


L pa = 1 d f* f(x) 
"D0 ra—aal Goor™ 
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1.1.3 Caputo derivative 


Since Riemann-Liouville fractional derivatives failed in the description and modeling of some 
complex phenomena, Caputo derivative was introduced in 1967. The Caputo derivative of 


fractional order a (n— 1 <a <n) of function f(t) defined as 


Cpa ee : DET) - 
DO = nama f, Paap 


if 0 <a <1 then the Caputo derivative of fractional order a defined as: 


Garlens .._ 4 Ya @ . 
DOH) = rama f GaP 


1.1.4 A new fractional derivative of Caputo (The inverse of the new 


Fractional Derivative without Singular Kernel) 


The authors of the article [6] suggest a new definition of fractional derivative, which assumes two 
different representations for the temporal and spatial variable. The first representation works on 
time variables, where the real powers appearing in the solutions of the usual fractional derivative 
will turn into integer powers, with some simplifications in the formulate and computations. In 
this framework, it is suitable to use the Laplace transform. The second representation is related 
to the spatial variables, thus for this non-local fractional derivative it is more convenient to work 
with the Fourier transform. The interest for this new approach is due to the necessity of using a 
model describing the behavior of classical viscoelastic materials, thermal media, electromagnetic 
systems, etc. In fact, the original definition of fractional derivative appears to be particularly 
convenient for those mechanical phenomena, related with plasticity, fatigue, damage and with 
electromagnetic hysteresis. When these effects are not present it seems more appropriate to 
use the new fractional derivative. We have also proposed a new non-local fractional derivative 
able to describe material heterogeneities and structures with different scales, which cannot be 
well described by classical local theories. So that, we rely that this spatial fractional derivative 
can play a meaningful role in the study of the macroscopic behaviors of some materials, related 
with non-local interactions, which are prevalent in determining the properties of the material. 
This work also contains some applications and simulations related to the behavior of these new 
derivatives, applied to classical functions such as trigonometric functions. These simulations 


show some similarities with the corresponding results by usual fractional derivative. 


16 


CHAPTER 1. FRACTIONAL CALCULUS 


Let us recall the usual Caputo fractional time derivative (UFD, ) of order a, given by: 


Oren 1 iG : 
DOH) = Tay f GaP 


with a € [0,1] and a € [—co,t), f € H'(a,b), and b > a. By changing the kernel (¢ — 7)~° with 


the function exp -—1 and Ti —a) with i 


They obtained the following new definition of fractional time derivative NFD; 


Definition 1.2. For a € [0,1], the fractional time derivative 9™u(x) of order (w) is defined 


as follow: 
Q f(x) = me) / f'(s) exp “=| ds (1.1) 
LF fe l-a 
where a € |— 00,2), b >a, f € H'(a,b), and M(a) is a normalization function such that 


M(0) = M(1) = 1, for example M(a) = 1 —0.2sin(27a), as shown in fig. 1.2. 


Figure 1.2: An example of a normalization function. 


The following Matlab function M has used in our programs: 


1 function y=Ma(z) 


2 r=0.2; 
A y=l4trx«sin (2* pixz) ; 
4 end 


The NFD, is zero when f(t) is constant, as in the UFD;, but contrary to the UPD;, the 
kernel does not have singularity for t = 7. Formulated also NFD; can also be applied to 


functions that do not belong to H'(a,b). Indeed, its can be formulated also for f € L'(—o00, b) 


ie 
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and for any a € [0,1] as 


at — aM (a) [ _ alt —T) 
AF) =F fH) - FO) en |-F— a 
Now, it is worth to observe that if we put 
l-a 1 
Oo= A €[0,c0o],a= 77 € 1G] 


the definition of NFD; assumes the form 


GO p(t) = = @) / “f(r)exp -— dr 


this is definition on which our study is based. 


Definition 1.3. [fn > 1, anda € [0,1] the fractional time derivative ger”) F(t) of order 
(n+ a) is defined by 
Be) $(t) = DE (BP FO). 


Theorem 1.1. For NF'D,, if the function f(t) is such that 


f(a) =0,s =1,2,..,n 


then, we have 


A” (FO) = 9 (BPO) 


Proof. [6]. We begin considering n = 1, then from definition (2.8) of ger F(t), we obtain 


Ro (Ar) = ue) [ f(r) exp -) dr 


1 


Hence, after an integration by parts and assuming f’(a) = 0, we have 


A" (910) = [ (Zr) ew-2ar = 


na | [ 7 (/ @) exp tar 


a a a(t — T) 
ert f'(r) exp— i ir 


a f@- — / f(r) exp ar 
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otherwise 


dt \l-a —a 
“To ro- fe pat ar] 


It is easy to generalize the proof for any n > 1. 


Definition 1.4. Let n > 1, anda € [0,1] the fractional derivative 9°*" f of order (n+a) is 


defined by 
Be F(t) = BE (9 FQ) 
= ue) [re exp - 2) ds 
Such that 
(a+n) M(q) : (n+1) a(t 7 s) 
RB joys = hi (s) exp | —_ | ds (1.2) 


The following MATLAB code calculates the new Caputo derivative for any C"*!-continuous 


function. 


1% The goal of this program is calculate 


2% the new Caputo derivative 
3 cle; close; clear all; 
4 syms t Ss; 


5 a=0; b=7; 

6 alpha=1.5; 

7 P(t =" 1254 

8 K(t )=Da(f ,alpha ,a) 


1 function y=Da(f,alpha ,a) 


9 syms t s; n=floor (alpha) ; 

3 y=int (diff(f(s) ,m+1)*exp((alpha-m) *(s-t) /(m+1-alpha) ) 
4 ,»s,a,t)*Ma(alpha -m) /(m+1-alpha) ; 

5 end 
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Lemma 1.1. Let n > 1, andy € (n,n +1). If we assume u € C"*'la, bl, then the fractional 


deferential equation 


Q” u(x) = 0, Vx € [a, b] 


has u(x) = Ant" + An12" 1} +++» + a,x +0; a; € R,i =0,1,--- ,n, as unique solutions. 


Proof. 0. Let u(x) = anx™ + Gna" 14+ +++ a,2 + a9; a; € R,i = 0,1,--- ,n, we have 


ul) (¢) = 0 such that Y™ u(x) = 0, Va € [a,b]. 


0. Let y €]n,n+ 1], it can be written of the form: y = +a where a €]0, 1[, and n = [4], 


we suppose J u(x) = 0, Vx € [a,b], we have 


Girly Oe) | * ail) (5) exp -e*) ds 


~ L=a 


and the Leibniz integral rule gives the formula 


rma @))= ute) Ce) 
Map” ess a(z — s) 
| ee p|- a J as 
= A) men) — 29 u(2, 


and we have J u(x) = 0, consequently ut (t) = 0 therefor, 


u(x) = Ant” + Qniz” 1 +---+ay2 +49; a; € R,i=0,1,--- ,n. 


The proof is complete. 


Lemma 1.2. Let y € (n,n+1),n = [y| > 0. Assume that f € C"\a, bl, the fractional derivative 
of order (y) belongs to Cla, bj. 


Proof. 0. Let a € (a,b) and s € (a, 2], we use part integration to prove this relation 


Simply to say that a function Y™ f is continuous (a, b). 
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0. Then, let {x,},,c, real sequence define in [a, b], which this sequence converges to x € |a, )). 


— MO) gm ( ies on) 


l—-a as 
(ay lim In —@ 
— f(a) exp | - a a ! 
lim x a( lim z,—s 
Mere im &n (n) (tim, u ) 
= FRG) free” $0 (e)exp |- | as, 
a-M(a r (n) came) 
ceri es| ae | . 
= 9 (a) 


Therefor, Y™ f is continuous [a, b]. 


1.2 Differentiation under the Integral Sign 


Suppose a function ¢ is given by the formula 


d 
ala) = J Flo, hat eis ea 


where c and d are constants. If the integration can be performed explicitly, then ¢’(x) can be 
found by a computation. However, even when the evaluation of the integral is impossible, it 
sometimes happens that @’(a2) can be found. The basic formula is given in the next theorem, 


known as Leibniz’ Rule. 


Theorem 1.2. Suppose that @ is defined by 


d 
ale) =f Fle, that es | 
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where c and d are constants. If fand f, are continuous in the rectangle 


then 
d 
¢' (x) = / ff tat, a22<b 


That is, the derivative may be found by differentiating under the integral sign. 


We consider a function defined by 


u1(2) 
(x) = / f(,t) dt (1.3) 


where uo(xz) and u;(x) are continuously differentiable functions for a < x < b Furthermore, the 
ranges of up and u; are assumed to lie between c and d. To obtain a formula for the derivative 
@'(x), where ¢ is given by an integral such as eq. (1.3), it is simpler to consider a new integral 


which is more general than eq. (1.3). We define 


F(g4),2) = [ 7 eta, (1.4) 


and obtain the following corollary of Leibniz’ Rule. 


Theorem 1.3. Suppose that f satisfies the conditions of Theorem (1.2) and that F is defined 
by (1.4) withe < y, z<d. Then 


a : i Alun vae (1.5) 
> Ge) (1.6) 
OF 

ae =f (a2) (1.7) 


Proof. Formula (1.5) is Theorem (1.2). Formulas (1.6) and (1.7) are precisely the Fundamental 


Theorem of Calculus, since taking the partial derivative of F' with respect to one variable, say 


y, implies that x and z are kept fixed. 


Theorem 1.4 (General Rule for Differentiation under the Integral Sign). Suppose that f and 


Of /Ox are continuous in the rectangle 
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and suppose that uo(x),ui(x) are continuously differentiable for a < x < b, with the range of 


Uo and u, in (c,d). If d is given by 


u1(2) 
1) = 4) d 
o(c) i f(w,t) at 
then 
u1(2) 
d(x) = f |x, m1 (2)] u(x) = f{2,uo(a)] -u(0) + f f(x, t) dt. (1.8) 
uo (x) 


Proof. We observe that 
F (x, uo(#), u(x)) = (2), 


in Theorem 1.5. Applying the Chain Rule, we get 


(2) = Fr + Fyup(x) + Feu} (2) 


Inserting the values of F),, Fy, and F, from eqs. (1.5) to (1.7), we obtain the desired result 
eq. (1.8). 


Example 1.1. Find ¢/(x), given that 


2 


- t 
Ox) = | arctan —dt 
0 a 


Solution. We have 
: —2t/x° 2tx 
— { arctan = =— 
Ox a 1+ (¢?/2*) eae 


We use formula (1.8) and find 


* Otxdt 
t? + x4 


¢'(x) = (arctan 1) - (2x) — [ 


Setting ¢ = xu in the integral on the right, we obtain 


ge) = f eer me, fl 2udu = 2 (2 —Iog2) 
0 0 


2 cu? + x4 , 


Theorem 1.5 (The Cauchy Formula For Repeated Integration). The Cauchy formula for 
repeated integration, named after Augustin Louis Cauchy, allows one to compress n anti- 


differentiation of a function into a single integral. Let f be a continuous function on the 
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real line. Then the n-th repeated integral of f based at a, 


=f [- fs (On) don +++ dogdoy 


is given by single integration 


Proof. A proof is given by induction. since f is continuous, the base case follows from the 


Fundamental theorem of calculus 


where 


Now, suppose this is true for n, and let us prove it for n +1. Apply the induction hypothesis 


and switching the order of integration, 


ey a= f I of Pil Ont \donsi-: -dagdo, 
a oa sf ce (o, —t)" f(t)dtdoy 


-— f (01 — 1)" f(Q\doyat 


=+ fw —ayrpqat 


n! 
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In this chapter, we interest concerns an integro-differential equation of Volterra when the 
derivative of the solution appears under the sign of integration in a non-linear way. The 
derivation is conformable fractional of Caputo type, which increases the interest of this study. 


We built Lipchitz conditions to obtain the existence and uniqueness of the solution in 
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a Banach space specific to the derivative of Caputo. The analytical study is followed by a 
complete numerical study. 

It should be noted that most of the integro-differential equations which have already been 
studied in case of the fractional derivative of the unknown u is out integral sign as shown in the 
references [6, 21], or in case the derivative of the first degree of the unknown under the integral 


like in the articles [45, 17, 46]. The authors of [45] studied the problem from the following type: 
wr) =F) +f K (a, s,u(s),u'(s)) ds, Va € [a, 6). 


where, f € C!({a,6]) and u is unknown, to be found in the same space, and in the paper [17], 
the authors studied for an integro-differential nonlinear volterra equation with weakly singular 


kernel as fellow 


all) = / p(t — s)K (t, s,u(s), u’(s)) ds + f(t), Vt € [a, b].. 


Therefore, the object of this chapter is to study a nonlinear integro-differential equation of 


Volterra of the type: 


f(x) = g(x) + [« (23,7 (8), QZ f(s)) ds, Va € {a,b}. (2.1) 


Where J f is the new of fractional derivative of Caputo (NFD), without singular kernel (see 
[6]), and this derivative has many properties as in the reference [27]. 

The first work, we built Lipchitz conditions to obtain the existence and uniqueness of the 
solution in a Banach space specific to the derivative of Caputo, then we search a method to 
approximate this solution and estimate error. In the last part, we give three numerical examples 
to illustrate the above methods for solve this type integro-differential equations. The exact 
solution is known and used to show that the numerical solution obtained with our methods is 


correct. We used MATLAB and C++ to solve these examples. 


Definition 2.1. We consider the space of continuously differentiable functions 


FQ (a,b) ={f €C(a,b), || f lla< co}, (2.2) 
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FQ“ (a,b) is a Banach space specific to the new Caputo derivative admits the norm 
| Flo = mane [f(a)|-+ max |B F(2)] Il f loo +I BF lo (2.3) 


Proposition 2.1. For a € (0,1), if f is continuous on [a,b], then the function h(x) = 9 f(x) 


is continuous on the interval |a, b]. 


Proof. 0. Let a € (a,b) and s € (a, 2], we have that 


# (sore [-2]) = Pore [2] +72 rere [-22=9], 


1 l—a Ila 
and hence 

poenn[-8E=B) - $ (se) e9[- AEB) ~ 2 rede [-2E=9) 
Then 


2 
sR 
So 
lo) 

oa 

a 
oa 
ell 
8 

| 
ey 
et 
Q 
H 

II 
a; 
| & 
—~ 
= 
oF 
lo) 
| rs | 
ele 
| 8 
2| « 
—% 
NY 
Q 


as a result 


oe DOT iets -“e= 2) ds. (2.4) 


Simply to say that a function J f is continuous (a, b). 
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0. Then, let {,},,c, real sequence define in [a, b], which this sequence converges to x € {a, )). 


lim 9 flan) = im ue) Ga, “exp -—"))) 
n—+0o —-a l-a 
a-M(a (Ln — 8) 
lm Toae [foe ts 1 | ie 
lim t,—-a 
M(a) is (lim : ) 
~ Ta f (im a) ~ F(a) a 1—a 
lim « a| lim z,-—s 
M (a) a(x — a) 
A (s(0) — Fa) exp |- 
a-M(a) [* a(xz — s) 
= 9) f(z). 
finally © f is continuous on [a, b]. 
Proposition 2.2. FY™ (a,b) is a Banach space. 
Proof. It is obvious that || - ||, is a semi-norm, and since ||- ||. is strictly positive, and || - || > 


ll - lo, we see that || - ||q is strictly positive too. Now, suppose that (fn),cn C FI (a,b) isa 


Cauchy sequence. 


Then, 
fn — fll, 79 asn,m— oo 
and hence 
fn — fmlloo 79, |B fr-I™ fnl|,, +0, as n,m — co 
Thus, the sequences (fn), and (9 eee are Cauchy sequences in 
(C({a, 6], R), || - ||.0) Since this space is complete, both sequences have limits in this space, we 


name these limits f and h, i.e. 
lfn — fll, 7 0 and a fn— All, 40, as n+ 0. 


The proof will be complete if we can show that f ¢ FI™ (a,b) 
and lim ||f, — fll, = 0. 
n+ +00 
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From eq. (2.4) we have 


l-a 


B f(0)+ F(a) -exp |-“E=9) 4 = [ 10)e |-"F =) ae. 


we note that for every n € N, 


Then 


lf - FL = sup |fa(@) - F@)| 
oo x€ [a,b] 


to su (a) tv) x 
< aay ow, | fa(x) — h(x)| 
acs, Pepto erate | Se) 
a sf ees _ a(t ~ s) 5 
“malts [0-10-42 


<C1||B fa — hlloo + Call fn — flloo + Callfn — flleo —> 0 


’ 
n—- +00 
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and 


9° fA], = sup, [9 Ale) — Ae) 
< TS sup Wale) — Fe) 
ey it -oya[-2= 
sap | [ tns)- fee || as 


< Cillfn 7 Filles + Callin an F lkee + Galli = Files ae 0, 


i.e., fy, converges to f and Y™ f,, converges to h in the uniform norm. Since the limit of a 


sequence in normed space is unique, we deduce that f = f and h =h. Then 


Ilfn — flla = Win — floc + || I fn -— 9 FI I 
= Ilfn — fille + |a i hl. as. 


n—-+00 


Finally FY“ (a,b) is a Banach space. 


2.1 Hypotheses 


We consider the following nonlinear integro-differential equation of Volterra: 


f(x) = g(x) + [« (x, s, f(s), Qo f(s)) ds, Va € {a,b}. (2.5) 


where f € FY™ (a,b) and f is unknown, and the function 


K: (ab? R? 3 R 


(z,s, u,v) > K(az,s,u,v). 


is satisfy the following assumptions: 


a) — €C([a, }]*R?) 


_ K 
bam >0,V 2,5 € [a,b], Vu,v € R, max (|Foesue) 


se (a, s.u0)1) <7 
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eq. (2.5) bearing several information for the solution f. If both sides of the equation are 


derived (using Leibniz’s rule for differentiation under the integral sign) we obtain Vx € [a, }): 


l-a 


BQ f(x) = YO g(x) + Mla) [ < [« (s, f(r), 9 F(r)) ar] exp -¢ — =) ds, 


= 9° ga) + FD PK (s,s, 410),9 He) exp |-“E—"] as 


l-a 


+ [ Fon IO), 9% f(r) arexp -“e= 2) = 


l-a 


= 9°) gx) + TD fC (5,5, £66), 9° H(8)) ep [-—)]) as 
+ Ma) [ Gi “ (s,7, f(r), 9 f(r)) ar) exp -e= 2) ds, 
So that 
B fe) = 9 ga) + ZO) I (i (5,5, £66), 9" 1(8)) ep [-“=")]) as 


ie ae) [ ([ “ (s,7, f(r), 9™ f(r) ar) exp -2e — =) ds, (2.6) 


finally, we obtain 


Definition 2.2. For g ¢ FY™ (a,b), we define the functional W, by: 


U,: FY (a,b) — FY (a,b) 


Y~ — W,(¢) 
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where 


W(y): [a,b] —>R 
b> Wyly)(a) = gl) +f K (0,5,9(s), 9 9(s)) a 
Proposition 2.3. For g € FY (a,b), the functional VW, is continuous from FY (a,b) to 
itself. 


Proof. Let g € FY (a,b). it is obvious that V,(y) is continuous on (a, b). 
For all x € (a,b), 


F WA(g)(e) = 9 gla)-+ PP [ (K (5,5,0(6), 2 ols)) ex |-“==") ) as 


_ Me) [ ([ “a (s,7, 0(r), 9 y(r)) ar) exp -“* ds, 


l-a 


JY V,(y) is a continuous function on (a,b) and bounded by the constant ||g||q + ABm(b — 
a)(b—a+1), where 


\= max we) and 6= max (exw -“=*)). 


a€(0,1) lL -—a@ ax<s<a<b l-a 


therefore U,(y) € FA™ (a,b). 
Let {Yn},,cn be a sequence of real-valued functions of the space F J“ (a,b) which converge 


to a function y € FY™ (a,b), 


Kim _Wy(n)(z) = tim (fH (ts. (s), 2 en(s)) ds + (2) ) 


n—++00 n—+-+00 


x 


— lim K (x, 8, Yn(s), 9B Yn(s)) ds + g(x) 


n—-+00 
a 


= [ lim (x, 8, n(s), 9 ¥n(s)) ds + g(x) 


n—-+00 


-/ K (145, lim p(s), lim I gals) ds + g(x) 
5 N— +00 nN— +00 


: [ K (2,8, 9(s), 9 p(s) ds + g(x) 


= V,(¢)(2). 


In this article, the purpose is to get the existence and uniqueness of the solution of eq. (2.1), 
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but use the minimum conditions to ensure its. 


2.2 Analysis study 
Theorem 2.1 (Existence). The equation (2.1) has a solution in FY™ (a,b). 


Proof. We define the following set: 


B= yp FI (a,b),| |Wy(g)(a) — f(a)| < (b- a)m, 


|IO U,(y)(xz) — BO f| < ABm(b— a)(b- a +1) 


It is obvious that, F is convex and closed set. For all y € E and all x € [a,b]; V,4(y)(a) = 
g(a), and 


No(eV2) - 92)] =] K (5,019), v(s)) as 


l—a 
He PS @ rs. 9 10) =| ds 


< \Bm(b — a)(b—a +1). 


Then, U,(E) C E. But, WV, is continuous on FY“ (a,b) to itself, and for all x, x’ € [a, }] 


(a <a") 
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We use the theorem Schauder for we proved that the functional VU, has a fixed point in 


FQ“ (a,b). Which result that eq. (2.1) has solutions in FY“ (a,b). 


2.3. Uniqueness 


Lemma 2.1 (see [45]). [f the function h is continuous and positive on (a,b), and satisfies 


an > 0, hie) < mf h(s) ds. (2.8) 


Then Vax € [a,b], h(x) = 0. 


Proof: (see [45]). Since h(a) is continuous function on [a,b], 4 « > 0 realizes the following: 
h(x) < p, Vt € [a, b]. 


Then, 


we integer 
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If we repeat this operation n times, we get: 


In order to prove the uniqueness of solution of eq. (2.1), we use the hypothesis: 


(c) 4My, Mo, M1, Mz € Ry, Vu, v,U,0 € R, V2, s € [0, UJ, 


(Hz) |K (x, 8, u,v) — K(z, 8,4, v)| < Mi|u —u| + Malu — 9 


OK OK —_ _ 
=— (2,5, 0,9) < M,lu — 4 + Mov —v| 


Theorem 2.2 (Uniqueness). The solution of eq. (2.1) is unique. 


Proof. 
Let f,h € FY™ (a,b) two different solutions of eq. (2.1). we have: 


and 


l-a 
= (Sora 108) 92 


: l—a 
+ we) / ( : a“ (s,7, h(r), BD h(r)) ar) exp -“= 2) ae 
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we define 


R(x) = |f(x) — A(x)| + |9© f(x) — 9 h(a)| 


oi< f IK (c, 5, f(s), 9 f(s)) — K (e, 8, h(s), 9 h(s))| ds 
i, (Mi|f(2) — h(2)| + Mo |B f(x) — 9 h(2)]) ds 
= 7 max(M,, M2) (|f(x) — h(z)| + |9© f(x) — 9 h(x)|) ds 
< / ” pani AB as 


< C1 | R(s) ds. 


we use the hypothesis (H2) 


A F(a) — 9 H(a)| <P] (5,5, 706). 9 £0) 
—K as s,h(s), 9 h(s)) | exp ia d 
$e) PEE (5.7, Fr), 9 F(7)) 
_ Ron. 1 f(r), B $(7))) ar-exp |-“E—" as 


<A im ine Ni Me Reet Ge) if max(M,, Mo) R(s) ds 


<C2- / R(s) ds. 


Consequently 
Gans max(Cy,C2) f R(s) ds 
Using the lemma 2.1 
lemma 2.1, 


R(x) <m [Rls \de— > Rig) =U; ve e la.0| 


we get: f(x) = h(x), Vz € [a, Bl. 
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2.4 Numerical study 


Under the assumptions (H,) and (H»2), we have prove that eq. (2.1) has a unique solution 


FJ™ (a,b), we introduce algorithm for finding numerical solution of eq. (2.1). Here, the 


interval [a,b] is divided in to N-equal sub-intervals, where x; = a+ j-l for all 7 € {0---N} 


and | = (b—a)/N. Also, using the following numerical integration formula: 


F: f(s)ds = UD wif (23) 


where w,; are the weights of the quadrature rule such that JW > 0, max, |w;| < W, on 
SIS 


eqs. (2.1) and (2.7) we get the system: 


fo = g(a) (2.9) 
ho = 9 g(a) =0 (2.10) 
hn = BO g(atn) + Tent owk (15,05, fj, hj) exp [-a 
+ reel tlt - a) (xj, 03, fj, hj) exp [-Sa (2.11) 
fn = 9(2n) 1D HK ny fi,hj). (2.12) 
ia 


Where f, approaches f(z,), and hy approaches J f (xp). 


2.5 System study 


In order to prove the existence and uniqueness of the solution of the system eqs. (2.9) to (2.12), 


we use hypotheses (H1) and (H2): 


Theorem 2.3. For | is sufficiently small, the system eqs. (2.9) to (2.12) has a unique solution. 


Proof. We supposing that the space 


R? has the first norm: 


U 
ER’, = |U|+1V| 


V 
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Now, we define: 


U Pl 
Yn > 1 Vv, == ; 
V P2 
where 
n-1 
pi = G(2n) +1 WK (an, Ln, U,V) +1) > wj)K (Pa tetas 
j=0 
M 
po = D™ g(an) + ; (2) W,K (Xn, Ln,U,V) 
—a 
M(a),< OK a 
+ To (« (5,25, fy hg) + (@n — 2A) Z (p25, fi hi) exp - 
joa 
We have 
U U' py 
Vn —_ Vn, — 
V y’ Bo 
1 1 
as result 
By = W, Uk (ase U, V) —kK (ius In; U’, V’)) 
M 
pia aa 12), We EG tat V = Git UG) 
however 


|21| <1 W (M,|U — U'| + M2|V — Vl) 
|G2| < Al W (M,|U — U'| + MQ|V — V") 
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Then 


<1 W(M,|U —U'| + MQ|V -V")) 
+ Al W (My|U — U"| + Ma|V —V’)) 


< max (WMI, W Mol,W M,Al, W M2Al) = 
V Vv! 


1 


Using the theorem of Banach to get a unique solution of the previous system eqs. (2.9) to (2.12) 


2.6 Error analysis 


In this part, we want to show that the numerical method constructed in the previous section, 


converges to the exact solution of the eq. (2.1), we define: 
€n = |fn — f(2n)| + Pin —~g9 F (20) 


We say that the method is convergent if 


lim ( max (en}) = () 
h—0 \ 0<n<N 
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For n > 0 and y € FY“ (a,b), we define the local consistency error by: 


O(h, Ln) a / K (i, Ss, y(s), go y(s)) ds = LS > w,K (Dail Oey h;) 


M(a) [ire s, (8), 9 (s)) exp [9 ds 


1D wk (25,05, fj, hg) exp [-s=)| 
j=0 


M (qa) 
l-a 


+ 


[e- _ (2.5,4(9).2 o(9)) en [-2F=9)] as 


l-a 


jee 
TL le j Be EG aed alt j) eXp - 


We say that the approximation method of the system eqs. (2.9) to (2.12) is consistent with 
eq. (2.1), if 
Vy € F9Y™ (a,b), lim (1 max <19(h, v»))) =) 
h—>0 \O0<n< 


Theorem 2.4. [f the approximation method eqs. (2.9) to (2.12) is consistent with eq. (2.1), 
then 


ale ( max, {60} ) = 
Proof. For n > 1, we have 
f (en) = gla») + i, K (0p, , f(s), 9 f(s) ds 
g(%n +1 wk (Daxtget 2h), OO F@,) 
and 
fr =G(@n) +1) wziK (tn 25, £5504) - 


So that 
If (@n) — fal = |E>_ w; (K (@n, 25, f(@;), 9 f(@;)) — K (Gn, 25, 5, hy) 
j=0 


< l ws, |i (tn, 05, f(x;), 9 f(2;)) > K Gy tinh) 
j=0 
<IWYS) (Mi |f(aj) — fy] + Mo |B F(a;) — Ay). (2.13) 


j=0 
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And 
(a) A(Ln — £;) 
hn Q Me eee Pe) Exp | — 
l-a 
a A(Ly — @; 
n— £j)—— i ti Fi j) €Xp - ( aad 


B® f(0) = 9 92) + 2) [i (5,5, 416). 9 t(9)) ex |-“E=| as 


M (qa) 
l-a 


[e ~ 3) (5,5, fs), 9 F(s)) exp [-2e=2) 


QB) f (an) = 9 g(a) + Ma) [- K (s,s, f(s), 9 f(s)) exp -= 4 ds 


+e Fen F (6.5, f0),9 M9) exp oa : 

~ YG g Ten owk (25,25, (23), 9 f(wj)) exp [-e 
M(a) a(Ln — £5) 
Tae l (Ln xi)a 25,56), 2 'F(«5)) exp | l-a | 
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Thus 
M ‘ ss 
|B flan) — bal & me), ae (25,07, F Gs), 9 f(s) exp <6 2) 
_ So wjK (7,23, fj, hy) exp [-a 


+ S/ wj(an — re (xj, 05, f(x;),9™ f(x;)) exp [-“a — = 


l-a 


~ OK 
— S70; (en - 25) (js Bir Fa, hej) EXP | 


j=l 


a(n 2] 


<M WSO (Mi |f(;) — fl + Mo [9 fay) — yl) 


j=0 
ube wy (M1, |f x3) — f)| + M2|9© f(@,) — hl) 
s 
finally 
|B fen) — hn| <M wih + (b—a)M)) |f(2;) — fj (2.14) 
= 
LW > +(Mz + (b— a)M2) |B f(x;) — hyl. (2.15) 
j=0 


From eqs. (2.13) and (2.14) we get: 


im, (amas. {en}) = 0 


h—0 \ 0<n<N 


2.7 Numerical result 


Since the approximation method of eqs. (2.9) to (2.12) converges to the solution of eqs. (2.1) 
and (2.7), we use the trapezoidal integration method to get that the terms f, and h,, will can 


not calculated exactly. but we use the iteration method of Banach for be approached its. 


Example 2.1. Let’s take the following nonlinear Volterra integro-differential nonlinear equation 


as an example: 


Va € [0,7]; f(x) = g(x) + [x (z,s, f(s), 9 f(s)) ds, (2.16) 
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where a = 0.5 and the kernel 


z+e*+sin?(s) + (e-* — V2sin(s + ay): 


ae ates ++ 


If we take: g(x) = sin(x) — x, we get f(x) = sin(z). 
The following MATLAB with C++ codes solve the equation (2.16) numerically, and plot 


the absolute errors and the approximate solution of this equation with the exact solution. 


1% The goal of this program is calculate the approximate 


2% solution of the following problem 
3 % u(t )=f(t)+\int(a,t)k(t,s,u,Du)ds t in [a,b] 
4 cle ; clear all; 


5 WIVTTIIIIIIN: Variables WQVIIIITIIIIVITI Do 
6 alpha=0.5; 

7 a=0; 

8 b=pi; 

9 PIVVVIVVTIVIVITN. functions SPYVVVVVTVTTVTSVIIITN 
10 syms t s x y; 

11 S(t)=txsin(t); 

12 DS(t )=Da(S, alpha ,a) 

13 k(t,s,x,y)= sin(t)*(t+texp(-s)+(s*sin(s))7*2 
14 +(exp(-s) - cos(s)+ sxcos(s) 

15 + sxsin(s))°*2) /(t+texp(-s)+x*2+y~°*2) 
16 K(t,s,x,y)= diff(k(t,s,x,y),t) 


1 #include ”moumen.h” 

2 double f(double) ,Ma(double) ,S(double) ,DS(double) , 
3 Df(double) ,KK(double ,double ,double ,double) ; 

4 double K(double ,double ,double,double) , 

5 k(double ,double ,double, double) ; 

6 main () 

7 { 


8 int N=64,n,j; 
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9 double alpha,a,b,t,h,beta,u[N+1], 

10 v [N+1] ,X[N+1],1[N+1],U,DU; 
11 alpha=0.5; 

12 a=0.; 

13 b=2.; 

14 h=(b-a) /N; 

15 beta=hxMa(alpha)/(1l-alpha) ; 

16 for (n=0O;n<N;n++) 

17 X[n]=a+hen; 

18 u[O|=f(a); 

19 v[O|=Df(a); 

20 1[0]=fabs(S(X[0])-u[0])+fabs (DS(X[0]) -v[0]) ; 
21 for (n=1;n<N;n++) 


22. { 

23 DU=KK(X[0] ,X[0] ,u[0] ,v[0]) *exp(alpha«*(X[0] -X[n]) 
24 /(1-alpha))/2.; 

25 Usk(X[n] ,X[0] ,u[0] ,v[0]) /2. ; 

26 for (j=1;j<n; j++) 

ay { 

28 DU-DUHKK(X[§] .X{4],u[5],v [i ]) exp (alpha (X[j] 
29 -X[n]) /(1-alpha)); 

30 UHK(X[n] ,X[5] uli] ,v(5]) 

31 } 

32 DU=DUHKK(X[n] ,X[n], f(X[n])+h*U, Df(X[n]) 
33 +betaxDU) /2.; 

34 U=U+ k(X[n] ,X[n], f(X[n])+h*U, Df(X[n]) 

35 +betaxDU) /2. ; 


36 v [n]=Df(X[n])+betaxDU; 
37 u[n]=f (X[n])+h«U; 
38 1 [n]=fabs (S(X[n]) -u[n])+fabs (DS(X[n]) -v[n]) ; 


41 FILE xfl ,*f2 ,* f3; 
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42 fl1=fopen(” Approximate.txt” ,”w”) ; 

43 for (n=0;n<N;n-+-+) 

44 fprintf (fl ,”"%lf \t Mlf \n*; meh 4 wl) )s 
45 fclose(f1); 

46 plot (” plot 'Approximate.txt' pointtype 5, 
AT x*xtx linewidth 2 ; 

48 pause mouse \n ”); 

49 endplot(); 

50 f2=fopen(” Error.txt” ,”w”); 

51 for (n=0;n<N;n+-+) 

52 fprintf(f2 ,”"%lf Vi Slt (n> mal . Lin) 33 
53 fclose(f2); 

54 plot1(” plot 'Error.txt' pointtype 7; 

55 pause mouse \n”); 


56 endplot1 () ; 


57 } 

58 double S(double t) 
59 { 

60 return txt-+t; 

61 } 

62 double DS(double t) 
63 { 

64 return 4.*t+2.*exp(-t)-2.; 
65 } 

66 double f(double t) 
67 { 

68 return txt; 

69 } 

70 double Df(double t) 
7m { 


72 return 4.*t+4.*exp(-t)-4.; 
73 +} 
74 double k(double t,double s,double x,double y) 
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79 


{ 


return (t+texp(-2.*s) -4.*s+4.*sxexp(-s)+8.*s*s 
+8.*s*s*xs + 4.*s*ses*xs + 1.)/(t 
4+2.*exp(-s)+4.*x*x+yxy/4.); 
} 
double K(double t,double s,double x,double y) 
{ 
return 1./(t+2.*exp(-s)+4.*x*x+tyxy/4) 
-(t + exp(-2.*s) -4.*s + 4.*sxexp(-s) 
+8.*s*s+8.*s*s*xs+4.%sxs*s*s+1. ) 
/((t+2.*exp(-s)+4.*x*x+yx*xy/4.) *(t 
+8. *exp(-s)+4.*x*xty*y/4.)); 


double KK(double t,double s,double x,double y) 


return k(t,s,x,y)+(t-s)#*K(t,s,x,y); 


double Ma(double t) 


{ 
double r=0.2; 


return 1l+rx*sin(2*pixt); 


} 
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‘Approximate.txt’ = 
sin(x) 


Figure 2.1: Numerical solution of example 2.1 for (N = 32) with the exact solution, and error 


obtained. 


N qaax, (lfi— f(ea)| + hi — 9 F(es))) 
32 0.01 

64 2.5E-3 

128 6E-4 

256 1.4E-4 

512 AE-5 

1024 8E-6 

2048 2E-6 

4096 1E-6 


Table 2.1: The numerical error obtained for the example 2.1 
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Example 2.2. Let’s take the following nonlinear Volterra integro-differential nonlinear equation 


as an example: 


Va € [0,7]; f(x) = g(x) + [x (x, s, f(s), 9 f(s)) ds, (2.17) 


where a = 0.5 and the kernel 


Qa + e—** + Asinh?(s) + 2(e* — 1)? 
nah ea Aue aa? 


Ke, 8,0, 0) = 


If we take: g(x) =e” —1—a, we get: f(x) =e” —1. 


60 


‘Approximate.txt' = 
exp(x)-1 


50 


40 


30 


20 


10 


Figure 2.2: Numerical solution of example 2.2 for (N = 32) with the exact solution, and error 


obtained. 


Example 2.3. Let’s take the following nonlinear Volterra integro-differential nonlinear equation 
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N max, (fi — f(vi)| + |hi — 9 F(#s)) 
32 L.8E-2 

64 4.5E-3 

128 1.2E-3 

256 3E-4 

512 7E-5 

1024 2E-5 

2048 4E-6 

4096 1E-6 


Table 2.2: Numerical error obtained for example 2.2 


as an example: 


—28 __ De — dg + Ase—* + 857 4+ 853 + 4st +142 


: ds, (2.18) 
x + 2e-* + (2f(s))? + (YO f(s) /2) 


Ve € [0,7]; f(z) e+] 
0 


where a = 0.5, we get: f(x) =t? +t. 


- gnax. (|fi— f(@)| + lh -— 9 F(2s))) 
32 4.5E-3 

64 1.2E-3 

128 1.2E-3 

256 2.5E-4 

512 ae 

1024 16 

2048 1E-6 


Table 2.3: Numerical error obtained for example 2.3 
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T 
‘Approximate.tt' 
X*X+X 


‘Error.ot' @ 


Figure 2.3: Numerical solution of example 2.3 for (N = 32) with the exact solution, and error 


obtained. 


In this chapter we prove the existence and uniqueness of nonlinear integro-differential equa- 
tion of Volterra with the use of the minimum of hypotheses that ensure this, and then we 
solved the problem numerically using numerical methods and programming the problem using 
the Matlab which is characterised by slowness, that needs the division of the program into two 
parts. A part using the Matlab and the second part using C++ which gave us satisfying results 


and in a short time. 
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A New Definition of Fractional Integral 
and its Properties with a Theoretical 


Application 


Summary 
3.1 The First Definition of New Fractional Integral ............ a9) 
3.2 Theoretical Application... ........... 0200000 ee eee 59 
3.3 The Second Definition of New Fractional Integral ........... 64 


In this chapter, we introduce a new definition of fractional integral as an inverse of the 
conformable fractional derivative of Caputo. By using these definitions, we obtain the basic 
properties of those fractional derivative and its fractional integral. Finally, we solve a class of 
fractional boundary value problems as a theoretical application, and we use Matlab to solve 
this class of fractional boundary value problems. 

The usual integral and derivative are (to say the least) a staple for the new technology, es- 
sential as a means of understanding and working with natural and artificial systems. Recently, 
many authors have participated in the development of the fractional calculus (differentiation 
and integration of arbitrary order). The applications of fractional calculus often appear in the 


fields such as generalized voltage dividers[50], electric conductance of biological systems|14], 
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capacitor theory|36], engineering|10], electrode-electrolyte interface models|27], feedback am- 
plifiers, medical[14, 39], fractional order models of neurons|4], analysis of special functions[24], 
and fitting experimental data|5]. 

In this chapter, we launch a new fractional Integral operator, we investigate some properties 
of the new fractional Integral operator. As concerns the properties of the fractional derivative 
operator, we are interested in recalling some extended functions 

As a theoretical application a class of fractional boundary value problems is solved. To 
the best of our knowledge, this is the first work that solve problem with the new concept of 
fractional derivative recently introduced by Caputo and Fabrizio in paper [6], which has many 


properties mentioned in article [27]. 


Lemma 3.1. Let a €|0,1[ and f € C'la,b] a non constant function. The equations 


Q™ u(x) = f(x), in [a,b] (3.1) 
wa) = FAA -ar'@) + aslo], in [a8 (3.2) 


have the same solution. 


Proof. If u is a solution of eq. (3.1), then 


MO) fv(syem [2] as = fe 


l-a l-a 


and the Leibniz integral rule gives the formula 


Ten paatea | Hew TS] =o 


MO cs 2 = ky 7 u'(s) exp |--¢ — a ds + f(z), 


eee ~ 1l-al-a l-a 
consequently 
M(a) a 
ul (e) = f(a) + Fe) 
So that 
ula) = Fay (Fla) + (1a) F@) 
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Example 3.1. Let a= 0.5, M(a) = 1, and f(x) = 12, the solution of the equation 


Q” u(x) =a 
is the same solution of the equation 
(@) = = (af (a) + (1a) f'(2)) 
ut (n) = Mla) af (x aly tz 
shut 
=5tt5 


Therefore 


1 1 
u(t) = 7a + sate, ceER. 


Lemma 3.2. Let y €]1,2] and f € C'la,b| a non constant function, then : The equations 


Q” u(x) = f(x), in [a,b] 


have the same solution. 


Proof. Let u be a solution of eq. (3.3), then, it satisfies : 


SE [ues - a ?) ne 


and the Leibniz integral rule gives the formula 


MO Dyna) —Q=Y MOD Pp) exp |- CVE“) ar = fe) 


2 ny 2-7 


M(y-1) , hee) y ! 
qa) = TB wl) + F'@) 


So that 
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Example 3.2. For y = 1.5 find the Green’s function for the following boundary value problem 


Q u(r) =2,0<2<1 


(3.0) 
u(0) =0,u(1) =0 
Solution: From Lemma 3.2 we have: 
1 
9” u(r) =a, 0< a <1eu"(s)= g(t + 1), O<2<1 
We get a boundary value problem equivalent to boundary value problem (3.5) 
" i 
u"(z) = ~(x + 1), O< a <1 
2 (3.6) 


The homogeneous equation wu” = 0 has the fundamental solutions u;(x) = x and u(x) = x—1 
to satisfy the boundary conditions Bo|u] = u(0) = 0 and B,{u] = u(1) = 0 respectively. 
Then W (ui, uz) (x) = 1 and therefore 


Thus solve eq. (3.6) with 


ua) =f sf(e)ds(a—1) + f (s —1)f(s)dsx 
=(e-1) fs-5(e+1)as+2 f (s-1)-5(s + 1)ds 


1 1 1 
_ A 223 Bo ee 2 
= -2 eo x x 2 a Ta 
= 1 3,14. 
oe Ge 37 
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Figure 3.1: Eexact solution of boundary value problem (3.5) 


3.1 The First Definition of New Fractional Integral 


In this section, we present our new definition as a theorem 


Theorem 3.1. Let n > 1, a € [0,1]. The formula : 


1 


IO = ara [ C-MoF9) + 1a) F(a 


where f € C[a,b], and M(a), is a normalization function such that 
M(0) = M(1) = 1 is a new fractional integral of order (n + a), and it’s as an inverse of the 


conformable fractional derivative of Caputo of order (n+ a). 


Proof. From Definitions 1.2 and 1.4, we obtain 


Bo f(t) = BO (9 F(t) 
ae / fOr) (s) exp |- — a ds, 


~ Tras —a 


using the following Leibniz integral rule 


b(t) 9 
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and considering a(t) = a, b(t) = t, and g(t, s) = f't!(s) exp | — ; 
—a 


= “ot =I, “ alt) = 0, and all s) = <a are) exp | we obtain 
Gas M 7 M en t— 
ten) Magog (Mf tag a) 
en, 
=Ylo+n) F(t) 
So that 
< (ge F(t) _ Mle) ‘aes (| = 7 ° — 9 f(t) 
We then obtain 
1 d 
ro) = ee [a a) 5 (BMI) Ha9!ef(0)]. 


Now, using Cauchy formula for evaluating the (n + 1)'” integration of the function f+ (t) 


f(t) ae [of aie (tn) dtn - ++ dtydt 
ee 


(¢—)"fO @dt 


nm! Jo 


= aan | (t—s)” c _ a) (9° f(s)) +aQet™ f(s)| ds, 
If we consider that g(x) = Be*” f(x) 
= 10 =F | - 9"lewls) + 1 a)g/(o)] ds 


Finally, we obtained an expression of the inverse of the conformable fractional derivative of 


Caputo as shown in the following figure: 


giatn) 
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The proof is complete. 


Lemma 3.3. Let y € (1,2). If we assume u € C'(0,1), then the fractional deferential equation 


Q” u(x) =0 


has u(x) = cx + C2; C1,Co € R as unique solutions. 


Proof. In order to prove the previous Lemma, it is sufficient to prove the following equivalence: 
Va € [a, 5], f(s)ds =0 Vz € [a, ), f(x) =0 (3.7) 
Suppose F’ is an anti-derivative of f, with f continuous on [a,b]. We have that 


F(x) — F(a) = [ f(s)as=0e¥e [a, b], F(z) = F(a), 


= Va € [a,b], F’(x) =0, 


= Ve € [a,)], f(x) = 90, 


Now, using (3.7) to prove Lemma 3.3, Vx € [0,1] we have that 


Gey =p 9 / WG ee | 7 a 7 a dr =0, 


is: Faces | wees a P=. 


= Va € [a, bj, u(x) = 0 


> Vax € [a, b], u(x) = crx +2; c1,c2 ER. 


Lemma 3.4. Assume that u € C1(0,1) with a fractional derivative of order y € (1,2), that 
belongs to C'(0,1). Then, those statements holds 


0. if u(a) =0, then 9M Fu(x) = u(z). 


0. IMM u(x) = u(x) +axt+b, a,bER. 


Proof. 0. Let y €|1,2[, we have that 


Lute) = 
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Q ula) = —— [ u' (rT) exp - © a a dr. 


Before outlining the method needed, we wish to recall the useful transformation formula 


[ [ [of te edie may f @- or see (3.8) 


For practical consideration 


/ : / " f(t)dtdt = / “(@ —1) f(t)dt (3.9) 


Using the formula (3.9) we obtain 


9” (Pu(ayy =I171 f fo =s)a(8) as dep - Sets a i 


~ 2-4 , dr 2-7 


+ PS [fe -sye06) 45] exp [-P PE) a 


7 7 u(r) exp - oi os = a d 
+ fue) exp - 0 oo a d 

_ — Pur) exp | O aes a: 2 d 
+ u(r) exp “ a = 2 

- = " u(r) exp O aa a d 
= u(x) — u(a) exp | 4 oa — * 


0. Firstly, we have that 


Q ula) = St [ ees - ileal a eo 
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We derive both sides of the equation, we get that 


< (9 u(x)) = we) 
7= 1 My = 1) i v= le=7) 
“5a, oe [w@en|- To9 Ja 
Mya) " y= il a7 
ia Sra (2) — 5 9 ula) 


using the last formula 


Py (g™ u(x) = UG=a) [« —s) C — Igy u(s) + (2— = (g™ us) ds 


3.2. Theoretical Application 


On the other hand, we will study the following fractional boundary value problem: 


(u(x) = f(x x 
9” ule) = f(z), 0S e<1 (3.10) 
“a(0) =O,a1)=0 


In the following, we present the Green’s function of fractional differential equation boundary 


value problem. 


Lemma 3.5. Consider a non constant function f € C(0,1), u is a solution of problem (3.10) 


if, and only if, it satisfies the following integral equation 
1 1 
in) = | Gi (x, s) f(s) ds +f Go(x, 5) f’(s) ds (3.11) 
0 0 
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where 
_ 4y-l s(a—1),0<a<s<l 
(y— 1) a(s—1),0<s<a2<1. 
and 
2 s(a—1),0<a<s<l 
Go(z, 8) = 7 


g(e— 1), On soe < 1. 


Proof. We may apply Lemma 3.4 to reduce first equation of problem (3.10) to an equivalent 
integral equation 


LQ” u(x) = f(x) + ax +b, 
for some a,b € R. Consequently, the general solution of problem (3.10) is 


1 


Way f ally DI) + @- YF W)lds+az+t 


u(x) = 
By second equation of problem (3.10), we get b = 0 and 


1 


a= ay f A -ally— DF) + @- FAs 


Therefore, the unique solution of problem (3.10) 


ule) = ayy ff @— llr — DAG) + @— Fas 
~ ay [0-9 - DF) + 2- VFA 
= ay ff W@- 9-20 ally - DFW) + @- Fas 
ey feb ilty- 46) + @-DFAas 
= amy fe Dilr- DFO) + @-V)FWlAs 
Tay [20 -oler- Die) +-aF'oHlas 


The proof is complete. 


Example 3.3. For y = 1.75 and M(a) = 1—0.2sin(27a), let find the Green’s function for the 
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following boundary value problem 


Q”v u(r) = 2? +2,0<2<1 
(3.12) 


Solution: From Lemma 3.5 we have: 


15 |}s(a@—1), 0O<a<s<l 


Gi(a, s) = — 
16 
a(s—1),0<s<a<l 
and 
5 |s(a—1),0<24<s<l 
Ga(z, $) = 76 


Thus solve eq. (3.12) with 


u(a) = f Gi(x, s)f(s)ds + Go(x, 5) f’(s) ds 


0 


15 =e 15 ft ' 
= (2 - = 2 
TAGs » | s(s ts)ds+ ie f (s — 1)(s° +s)ds 
5 . 5 f° 
+2@-1 f s(25-41)ds+ 2x f (s — 1)(2s + 1)ds 
16 : 16° J, 
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Figure 3.2: Eexact solution of boundary value problem (3.12) 


Example 3.4. For y = 1.75, let find the Green’s function for the following boundary value 


problem 


QV u(r) = 2? +2, 0<2<1 
(3.13) 


Solution: From Lemma 3.5 we have: 


15 |} s(a@—1), 0O<a<s<l 


Gi(z, 8) = 
16 
a(s—1),0<s<a<l 
and 
5 |s(a@—1),0<a4<s<l 
Ga(z, s) = 75 
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Thus solve eq. (3.13) with 


u(x) = f Gi(a,s)f(s)ds+ | Go(a,s)f'(s) ds 


0 


15 . bf 
ae aaa —] 2 paid —] 2 
TAGs ) f s(s ts)ds+ ie | (s — 1)(s* +s) ds 
5 on i 
+2@-1 f s(25-41)ds+ 2x f (s — 1)(2s + 1)ds 
16 ; a 


Finally, we program the problem (3.13) using Matlab to get the solution u(x), as follows: 


1% The goal of this program is calculate 
2% the new Caputo derivative 

3 cle ; close ;clear all; 

4 syms X Sj; 

5 a=0; b=1; gamma=1.75; 

6 f£(x)=x°2+1; 

7 I(x)=(int (s*(x-1) *((gamma-1)*f(s) 


8 +(2-gamma)*diff(f(s))),s,0,x)); 
9 J(x)=(int (x*(s-1) *((gamma-1)x«f(s) 
10 +(2-gamma)*diff(f(s))),s,x,1)); 


11 u(x)= expand(I(x)+J (x) ) /Ma(gamma- 1) 
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@ MATLAB R2019b - a x 
EDITOR VIEW s 
vo a LalrindFaes | Insert FE e al ~ SS > bal (&] Run Section i> 
oy Geen en Por Ceo ee S| a ae [Advance Runand 
. ¥ > G@Print ~ Q Find » Indent [E] $3] [Zé ~ ~ Advance Time 
FILE NAVIGATE EDIT BREAKPOINTS RUN z= 
S]> HF (| > D: > Doctorat » LaTeX3 +\2 
Current Folder @ Editor - D:\Doctorat\LaTex 3\myprotr.m Qa x 
Name « myprotr.m Ma.m [+ | 
« B G 
——— “1- cle ; close ;clear all; - 
{iil figl.PNG = 7 
| fig2.PNG 2 syms x s; 
a, 
ai ae ay= a=0; b=1;gamma=1.75; 
Brome, M4 — f(x) =x*24x; 
Gino ans = I (x) =(int (s* (x-1) * ( (gamma-1) *f (s) + (2-gamma) *diff(f(s))),s,0,x)); 
— Giz J (x) =(int (x* (s-1) * ( (gamma-1) *f (s) + (2-gamma) *diff(f(s))),s,x,1)); 
moumen - ART3.tex 
moumen-ARTAaux |) 77 — u(x)= expand (I (x) +J (x) ) /Ma (gamma-1) | 
moumen - ART4.log = 
“| moumen - ART4.out y 


Ma.m (Function) Command Window 


Workspace @ | New to MATLAB? See resources for Getting Started. 


u(x) = 


(5*x*4) /64 + (25*x*3)/96 + (5*x*2)/32 - (95*x) /192 


3.3. The Second Definition of New Fractional Integral 
In this section, we introduce a new definition of fractional integral as a theorem: 


Theorem 3.2. Letn > 1, a € [0,1], and f € C'la,b]. The formula : 


lft) = aa f (t= 9)" Malt 8) + mld — ad] f(s) as 


M (qa) 


where f € C[a,b], and M(a), is a normalization function such that 


M(0) = M(1) = 1 is a new fractional integral of order (n+ a). 


Proof. From Definitions 1.2 and 1.4, we obtain 


2 feel A 


l-a 


gour) f(x) = 


and Leibniz integral rule gives the formula 
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So that 
@ ¢glatn) ¢(4y) — M(® porta A glatn) 
ae) So ae a Oe 
We then obtain 
1 d 
ro) = ee [a a) 5 (Bem) Ha9ef(0)]. 


Now, we use the Cauchy formula for evaluating the (n+1)"" integration of the function f(*(t) 
; d 
10 = say [ ¢- 2" [a-9z (OOs@)) +0971) as 


Finally, we get 


F(t) = aay f 9M" lotls) + (1 a) Fo] es 
: eas —s)+n(1—a)|f(s) ds 
= Fa | 9 ats) + n(1 - ay) f(s) 4 


Lemma 3.6. Let y € (n,n+1),n = [7] > 0. Assume that u € C"/a, b], then those statements 
holds: 


0. Ty (9 u(t)) = u(t) + Dae Ge Re = Oise 
i=0 


0. if u(a) = 0, then QM (Iu(t)) = u(t). 


Proof. 0. Let y €]n,n + 1[, it can be written in the form: y = n+ a where a €]0, 1[, and 


n = [7], we have 


fi (g” ult) = ima | (t — s)" “a(t — s) + n(1—a)|9™ u(s) ds. 
= a (t — s)"B u(s) ds + af (t — ye (g™ u(s)) ds 
(1-a) y n{ M(Q) (nt a - 
* Ma) nl j ve) (a ae 'u() = 
1 t 


=— | (t—s)"u*(s) ds 


=u(t) + > ait’, a; €Ri=O0,1,...,n 
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0. Let 7 €]n,n + 1[, it can be written in the form: y = +a where a €]0, l[, and n = [4], 


we have 


B® (Kult) = 7 a i ae E / (e=a)'ulo)ae ae -“ | ds 
| | cas = i: ‘(e=a)'w (wae an - 2] ds 
2 aia f Hs es te] de 4 [uw ao |---| a 
C= 


= aia [ wen - | ds 
-ot- a9] - g [wren |-SE 9] ds 
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We study the existence and uniqueness of the solution of a fractional boundary value problem 
with conformable fractional derivation of the Caputo type, which increases the interest of this 
study. In order to study this problem we have introduced a new definition of fractional integral 
as an inverse of the conformable fractional derivative of Caputo, therefore, the proofs are based 


upon the reduction of the problem to a equivalent linear Volterra-Fredholm integral equations 
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of the second kind, and we have built the minimum conditions to obtain the existence and 
uniqueness of this solution. The analytical study is followed by a complete numerical study. 

Recently, papers have been published that deal with the existence and multiplicity of the 
solution of nonlinear initial fractional differential equation by the use of techniques of nonlinear 
analysis, see [36, 26, 39, 54]. However, most of the papers offer the problem using the stan- 
dard Riemann-Liouville differentiation. However, Our aim is to study the existence and the 
uniqueness of the solution for a class of fractional boundary value problems. To the best of our 
knowledge, this is the first work that solves problem with the conformable fractional derivative 
by Caputo and Fabrizio in paper [6], which has many properties mentioned in the article [27]. 
The interest for in this new approach is due to the necessity of using a model to describe the 
behavior of classical viscoelastic materials, electromagnetic systems, thermal media, etc. In 
fact, the original definition of Caputo’s fractional derivative appears to be particularly conve- 
nient for those mechanical phenomena, related to damage and with electromagnetic hysteresis, 
fatigue and plasticity. When these effects are not present it seems more appropriate to use the 
new fractional derivative [6]. 

In this chapter, We study the existence and uniqueness of the solution of the fractional 


differential equation boundary value problem, as follows: 


(4.1) 


where 1 < y < 2 is a real number, q is the potential function, and f : [0,1] — R is continuous. 
and J is the new fractional derivative, and we introduce a new definition of its fractional 
integral with some properties, using this fractional integral upon problem (4.1) to obtain an 
equivalent linear Volterra-Fredholm integral equations of second kind. Finally, by the means 
of some theorems, the existence and uniqueness of solutions are obtained, and we introduce an 


algorithm for finding a numerical solution of this problem class. 


4.1 Analytic Study 


In the following, we suppose the function M(qa) = 1. 
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Lemma 4.1. Given q € C[0, 1], and 1 < y < 2, the solution of 


Qo” u(x) + q(x)u(xz) = f(x), O<2x<l 


(4.2) 
u(0) = u(1) = 0 
satisfies the following linear Volterra-Fredholm integral equations of the second kind 
Ma 1 
u(x) +f G(x, s)u(s) ast f K (a, s)u(s) ds = g(x) (4.3) 
0 0 
1 


where g(x) = [i —1)(as —1+a)f(s)ds+ | xr(as — 1) f(s) ds, 


0 


G(a,s) =(x—1)(as —1+a)q(s) and K(x, s) = x(as — 1)q(s). 
Proof. We may apply Lemma 3.6 to reduce Eq. (4.2) to an equivalent integral equation 


Ij (9 u(z)) = Ip (F(z) — a(2)u(x)) 


>u(r)+cr+d= [iow —s)+(1-a)] (f(s) — q(s)u(s)) ds 


Using boundary conditions u(0) = u(1) = 0, we have d = 0, and 


= / (1 aa) (f(s) — ¢(e)ule)) de. 


Therefore, the unique solution of problem (4.2) is 


The proof is complete. 
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4.2 Existence and uniqueness of the solution 


The classical approach to proving the existence and uniqueness of the solution of (4.3) is the 


Picard method. This consists of the simple iteration for n = 1,2,... 
x 1 
t, (az) = 9(2) +f G(x, 5)Un_1(s) s+ [ K(a, 8)Un_1(s) ds (4.4) 
0 0 
with uo(a) = g(x). For ease of manipulation, it is convenient to introduce 
0,(£) =, (0) = tae), = 1,2)2- (4.5) 


with vp(a) = g(x). On subtracting from (4.4), the same equation with n replaced by n — 1, an 


we see that 
(2) = i K(@,4)0, (sds, a= 1,214. 
0 
Also, from (4.5) 


git) = dvi) (4.6) 


The following theorem uses this iteration to prove the existence and uniqueness of the solution 


under quite restrictive conditions, namely that G(z,s), K(2,s) and g(x) are continuous. 
Theorem 4.1. /f g(x) is continuous in0 < x < 1, and the function K(x,s), G(x,s) are 


continuous in0 <s <a < 1, and panax., A(x, s)| < 1, then the integral equation (4.3) 


possesses a unique continuous solution forO <a <1. 


Proof. Choose M,, Mz and M3 such that 


lg(z)|< Mi, O<a<l 


IG(z,s)|< Mo, O<s<aK<l 


|K(a,s)| << M3, 0<s<a2<1 where M3 < 1. 


We first prove by induction that 


noe 
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this bound makes it obvious that the sequence u,,(x) in (4.6) converges, and we can write 
ur) = > v;(x) (4.8) 
i=0 


We now show that this u(a) satisfies equation (4.3). The series (4.8) is uniformly convergent 
since the terms v;(x) are dominated by M,(Mox)'/i!+ M, M3. Consequently, we can interchange 


the order of integration and summation in the following expression to obtain 


[ G(z, s) S-ui(s) ds + ; Ke, &) S © ui(s) ds = ¥ [ Gle,s\u) ds 


i=0 i=0 i=0 
oo 


Each of the v;(x) is clearly continuous. Therefore u(x) is continuous, since it is the limit of 
a uniformly convergent sequence of continuous functions. 
To show that u(x) is the only continuous solution, suppose there exists another continuous 


solution u(a) of (4.3) Then 
u(x) — u(x) = [ G(x, s)(u(s) — u(s)) ds a K (a, s)(u(s) — u(s)) ds (4.9) 


since f(x) and f(x) are both continuous, there exists a constant C’ such that 
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For a large enough n, the right-hand side is arbitrarily small, therefore, we must have 


Theorem 4.2. /f f(x), q(x) are continuous in [0,1], and max lq(a)| < 1, then the fractional 


boundary value problem (4.1) possesses a unique continuous solution forO <a <1. 


Proof. If f(x), g(a) are continuous in [0,1], then it is clear that the following functions 


g(x) = [i —1)(as —1+a)f(s)ds +f x(as — 1) f(s) ds, 
G(a,s) = (x —1)(as —1+a)q(s), 


K(x, 8) = x(as — 1)q(s), 


are continuous, and |K (a, s)| = |x(as — 1)q(s)| < |g(s)| < 1, Vz,s € [0,1], which means that 


integral equation (4.3) possesses a unique continuous solution for 0 < x < 1. Therefore, there is 


a unique continuous solution of the fractional boundary value problem (4.1) for0 < x < 1. 


4.3 Numerical study 


In this section, we introduce an algorithm for finding a numerical solution of linear Volterra- 
Fredholm integral equations of the second kind, the methods based upon trapezoidal rule. For 
all N € N, Here the interval [0,1] in to N equal sub-intervals, where h = (b—a)/N, and 
x, =a+i-h for allie {0---N}. 


The formula of the numerical integration is: 


f(a) + 25> f(a) + FO) 


j=l 


f(s) ds = 


we apply this formula in eq. (4.3), and we obtain: 


G(a;,Xo)u(x%o) + 2 S Gly; )ule;) Gla, ata 


j=l 


g(x) = u(ay) + 


+ 


ml > 


K (xj, %o)u(a%o) + 2 ‘ I (is, 0; UG, ) = ae, ex)u(on) 


j=l 
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. h i-1 
=>Vi= 0,...,N, Gi = Ui t+ 3 pee eee + Giyu; 
h N-1 
es Kiguo + 2 d, Ayjuj + Ky 
This leads to 
h — hf 
3 (Gio + Kjo)uo + hS (Gi ae Crea 5 € PG 2K) Uy 
j=l 
N-1 h 
+h oa Ayjuj + 3 hintn = 9 
j=it+l 
Finally, we get a system of N + 1 equations, which is: 
AU =B (4.10) 


when B = (90, 91,---,9nw), U = (uo, t1,--., Un), and A = (Giz) ag ioe N3 


-Koo/2+1 ifi=j=0, 

Koj ifj=1,...,.N—-1, 

- Ko;/2 if 7 = N, 
(Chea)? ifi=1,...,N, 

(Gi +2K,)/24+1 ifi=j=1,...,N—-1, 


“(GgtKy)/2+l teagan, 


iS 
BS 
| 
ST TF TF F FF TF FF FTF FSF 


dons? ifi=1,...,.N—1, 
Ge ake) PESO NG SPs. PST 
ee eS NS ge Ny 
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K 1 K 
os Kou Koa es Kon-1 a 
Ki9+G K G 1 K. 
of 10 ut. He ; Ky Kin-1 4s. 
K20+G20 Koi+Go1 Ko2+Go2 a 4h : Kin 
2 2 2 h 7 2 
A=h 
Kn—-in-1+Gn_-1in-1 ns 1 Kn-1n 
2 h 2 
Kwnot+Gwo Knit+Gmi > 93 nes Knn-1+Gnn-1 Knn+Gnn , 1 
2 2 2 2 h 


We have chosen to write our system in its general matrix form without taking into account 


the fact that uo = uy = 0. However, we can see that 
go = Ko; = 0, Vj € 10.0} tig = 0. 


In same way, we get uy = 0. 


4.4 Numerical result 


In this section, we give three numerical examples to illustrate the above methods for solve the 
linear Volterra-Fredholm integral equations of the second kind. The exact solution is known 
and used to show that the numerical solution obtained with our methods is correct. We used 


MATLAB to solve these examples. 


Example 4.1. Consider the following Fractional Boundary value problem: 


Q”) u(x) + q(x)u(xz) = f(x), O<ax<l 
(4.11) 


where y = 1.5, q(x) =1, and 


_. SORE Be Se = ee 17x +17 
7 20 , 


f(@) 


with the exact solution u(x) = (e” — 1)(a — 1). 


The following MATLAB code solves the equation (4.11) numerically, and plot the absolute 


errors and the approximate solution of this equation with the exact solution. 
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1% ce programe qui calcul la solution approchee de la 


probleme 


OU) ake 


u(t)=f(t)+\int (a,t)k1(t,s,u(s))ds+\int (a,b)k2(t,s,u(s))¢ 


t in [a,b] 
3 cle ; clear; 
4 WGIYVTTIIIIIIS. Variables WY IISISISSSo 
5 gamma=1.5 ; 
6 alpha=gamma- floor (gamma) ; 
7 a=0; 
8 b=1; 
9 N=16; 
10 h=(b-a) /N; 
11 tic 
2 RVVTVVVVVVVIH functions WVVVVITVVTVTIIIN. 
13 syms t 8s ; 
14 u(t)=(t-1) *(exp(t) -1) ; 
15 q(t)=1+40xt ; 
16 f(t)=sinh(t) + (exp(t) - 1)*(t - 1) + txexp(t); 
17 V(t,s)=(alphax(t-s)+1-alpha)x«q(s) ; 
18 F(t,s)=t*(alphaxs-1)*q(s); 
19 g(t)=int (V(t,s)*f(s),s,0,t)+ int(F(t,s)*f(s),s,0,1); 


20 07070707070707070707070707070707070707070707070707070707070707070707070707070707070707070707070707070707070707070 


 0/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ OF O/ O/ O/ OF O/ O/ O/ O/ O/ O/ OF O/ O/ O/ O/ O/ OF OF O/ O/ O/ O/ O/ OF O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ OF O/ O/ O/ O/ O 


21 for n=1:N41 


22 X(n)=a+h«*(n-1) ; 
23 B(n)=g(X(n) ) 5 
24 end 


2 for i=1:N+1 


26 for j=1:N+1 

27 vi(i,j)=V(X(i) .X(j))3 
28 v2(i,j)=F(X(i) .X(j))3 
29 end 

30 end 
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A=zeros (N+1,N+1) ; 
for 1=1:N41 
for j=1:N+1 
v1(i,j)=V(X(i) .X(5)) § 
v2(i,j)=F(X(i) .X(j)) 5 
end 
end 
A(1,1)=1+h*v2(1,1) /2; 
ioe | =2:N 
A(1,j)=h*v2(1,j); 
end 
A(1,N+1)=h*v2 (1 ,N+1) /2; 
for i1=2:N 


A(i,1)=hx*(v1(i,1)+v2(i,1))/2; 
A(i,i)=1+h*(v1(i,i)+2*v2(i,i))/2; 


A(i ,N+1)=h*v2(i ,N+1) /2; 


end 


A(N+1,1)=h*(v1(N+1,1)+v2(N+4+1,1)) /2; 
A(N+1,N4+1)=1+h*(v1(i,i)+v2(i,i)) /2; 


for i=3:N+1 
for j=2:i-1 
A(i,j)=h«(v1(i,j)+v2(i,j)) 5 
end 

end 

for i1=2:N 
for j=i+1:N 
A(i,j)=hev2(i,j); 
end 

end 

U-vpa(inv (A) )*B'; 


for i1=1:N41 


r(i)=vpa(abs(U(i)-u(X(i)))); 
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64 end 

65 toc 

66 plot (X,U); 
67 grid on; 

68 plot (X,r); 


69 grid on; 


Figure 4.1: The Absolute Error of Test Example (4.1) with N = 16. 


Example 4.2. Consider the following Fractional Boundary value problem: 


QO u(x) + q(x)u(xz) = f(z), O<x<l 
(4.12) 


where y = 1.35, ga) = 1, and 


3(1697? + 49) sin(ax) + 78073(cos(ax) — e~ 7/13) — 4207? sin(r2) 


Ka) = 16972 + 49 


with the exact solution u(x) = 3sin(72). 
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Figure 4.2: The Absolute Error of Test Example (4.2) with N = 128. 


Example 4.3. Consider the following Fractional Boundary value problem: 


Q”) u(x) + q(x)u(xz) = f(x), O<2x<l 


(4.13) 
u(0) = u(1) = 0 
mere y= 1.75, g(t) = ee and 
20 
__ a(x + cos(mx))(2e~* + 19) 3 4n?(cos(1x) + e73* 
H(z) = 20 (TG Op 
i 8 _ e824 (32 + 827) cos(rx) + (x — 30)m sin(pix) + re~** 
3 (7? + 9)? 


with the exact solution u(x) = x(cos(7x) + 2). 
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-3 
x10 
3.5 


Figure 4.3: The Absolute Error of Test Example (4.3) with N = 16. 


Example 4.4. Consider the following Fractional Boundary value problem: 


(4.14) 


where y = 1.9, q(x) =t—1, and f(x) = 102. In this case, we don’t know the exact solution. 


The following MATLAB code solves the equation (4.14) numerically, and plot the approxi- 


mate solution of this equation. 


1% ce programe qui calcul la solution approchee de la 


probleme 
2 0. a8 
u(t )=f(t)+\int (a,t)k1(t,s,u(s))ds+\int (a,b)k2(t,s,u(s))4q 
t in [a,b] 
3 cle ; clear; 


4 WGYVITTIIIIIISo Variables WY TIIISIISSSISo 
5 gamma=1.9 ; 
6 alpha=gamma- floor (gamma) ; 


7 a=0; 
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8 b=1; 

9 N=16; 

10 h=(b-a) /N; 

ll WWVWTVIIITTH functions WTSI ITTIN. 

12 syms t s ; 

13 q(t)=t-1; 

14 f(t)=10«t ; 

15 V(t,s)=(alphax(t-s)+1-alpha)x*q(s); 

16 F(t,s)=t*(alphaxs-1)*q(s); 

17 g(t)=int (V(t,s)*f(s),s,0,t)+ int(F(t,s)*f(s),s,0,1); 


18 AV AVIV AVAV ADV AVAVAV AVAL A VAY AV AV AVAV AD DAV AVAV A AAVAVAV AV AVAVA VAD AVAVAV AA AVAVAV AV AVAVAVA DAV AVAVAAVAVAVAVAY 


0/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ OF O/ O/ O/ O/ O/ O/ OF O/ O/ O/ O/ O/ OF O/ O/ OF O/ O/ O/ O/ O/ O/ O/ O/ O/ O/ O07 O/ OF O/ O/ O/ O/ O/ O/ O/ O/ O/ O 


19 for n=1:N+1 


20 X(n)=a+h*(n-1) ; 
21 B(n)=g(X(n) ) ; 
22 end 


23 for i=1:N+1 


24 for j=1:N+1 

25 v1(i,j)=V(X(i) .X(j)); 
26 v2(i,j)=F(X(i) .X(j)); 
a7 end 

28 end 

29 


30 A=zeros (N+1,N+1); 
31 for i=1:N+1 


32 for j=1:N+1 

33 v1 (i,j )=V(X(i) .X(j))3 
34 v2(i,j)=F(X(i) .X(j)); 
35 end 

36 end 


37 A(1,1)=1+h*v2(1,1) /2; 
38 for j=2:N 
39 A(1,j)=h*v2(1,j); 


40 end 
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41 A(1,N+1)=h*v2(1,N+1)/2; 
42 for i=2:N 


43 A(i,1)=h*(v1(i,1)+w2(i,1))/2; 

44 A(i,i)=I1+h«(v1(i,i)+2*v2(i,i))/2; 
45 A(i ,N+1)=h*v2(i ,N+1)/2; 

46 end 


47 A(N+1,1)=h*(v1(N+1,1)+v2(N+1,1)) /2; 
48 A(N+1,N+1)=14+h*(v1(i,i)+v2(i,i))/2; 
49 for 1=—3:N+1 


50 for j=2:i-1 

51 A(i,j)=h*(v1l(i,j)+v2(i,j)); 
52 end 

53 end 


54 for i=2:N 


55 for j=i+1:N 

56 A(i,j)=h*v2(i,j); 
57 end 

58 end 


59 U-vpa(inv (A) )*B'; 
60 plot (X,U); 


61 grid on; 
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0.25 


0.2 


0.15 


0.1 


0.05 


Figure 4.4: The approximate solution of test Example (4.4) with N = 32. 
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Conclusion 


In this thesis, we first dealt with some Volterra integro-differential problems in fractional 
calculus by using the new derivative of Caputo, and we proved the existence and uniqueness of 
each solution of this problems. Convergences of the obtained solutions are also justified in order 
to establish that the formal solutions are analitic solutions. Afterwards, we have examined with 
a conformable fractional derivative of Caputo and its Fractional Integral. Considered problem 
reduced to the equivalent linear Volterra-Fredholm integral equations of the second kind. 

As perspectives, we will try to study the generalized Fractional Boundary Value Problem 
with with a conformable fractional derivative of Caputo and its Fractional Integral. This 
generalized version represents a challenge from the analytical point of view, i.e. the existence 


and uniqueness of the solution. However, the numerical side remains the same. 
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Annex 


MATLAB Operators and Special Characters 


Some common commands are listed in this section a full specification of each can be obtained 


using the help system. 


+ Addition. 

- Subtraction. 

= Multiplication. 

/ Division. 

abs Absolute value. 

sqrt Square root function. 

cle Clear command window. 

clear Clear variables and functions from memory. 

x Exponentiation. 

pi the mathematical constant 7. 

' Transpose. 

% To comment out one line in a multiline command. 
cle clean command window. After this function, all previous command written on 


window will clean. 
close all closes all figures, and window. 
clear all this build in function clear all variable created in work space of matlab. 
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syms Create symbolic variables and functions. 


diff(f) differentiates f with respect to x, and there are several forms 
diff(f(x)) = f'(x) 
diff(f(x)n) = f(z) 
; dd” x, 
diff(fy)my) = ew) 
y 
: d d 
diff(f,x1,...,.xN) = j= ge ices 
int Integrate, and there are several forms 


int(f,s,a,t) = [1 ds, int(f,omega) = | f(s) ds. 
a Q 


Da(f,alpha,a) the new fractional derivative of Caputo J f(t) 


floor the nearest integer in the direction of negative. 
Ma M the function used in section 1.1.4. 
A (i,j) A;; Element of the matrix A. 

V(i) V; Element of the vector V. 

== Equal to. 

v= Not equal to. 

> Greater than. 

>= Greater than or equal to. 

= Less than. 

— Less than or equal to. 

& Logical AND. 

— Logical OR. . 

n Logical NOT. 
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Compare Matlab to C/C++ 


Pros Matlab: 


e is a higher level of programming language. 

e is easier to start. 

e provides convenient tools and built-in functions. 
e provides user-friendly graphical interface. 

e is good for post data analysis and visualization. 
e Matlab is better for numerical treatments. 

e Matlab is more popular in engineering. 


e Symbolic Math Toolbox™ provides functions for solving, plotting, and manipulating sym- 
bolic math equations. You can create, run, and share symbolic math code using the 
MATLAB®) Live Editor. The toolbox provides functions in common mathematical areas 
such as calculus, linear algebra, algebraic and ordinary differential equations, equation 


simplification, and equation manipulation. 


Cons Matlab: 


e in many cases, the computational speed of Matlab is slower than that of C/C++. 
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